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Abstract 

Above-threshold ionization (ATI) results from strong field laser-matter interaction and it is 
one of the fundamental processes that may be used to extract electron structural and dynamical 
information about the atomic or molecular target. Moreover, it can also be used to characterize the 
laser field itself. Here, we develop an analytical description of ATI, which extends the theoretical 
Strong Field Approximation (SFA), for both the direct and re-scattering transition amplitudes 
in atoms. From a non-local, but separable potential, the bound-free dipole and the re-seattering 
transition matrix elements are analytically computed. In comparison with the standard approaches 
to the ATI process, our analytieal derivation of the re-scattering matrix elements allows us to study 
directly how the re-scattering process depends on the atomic target and laser pulse features - we can 
turn on and off contributions having different physieal origins or corresponding to different physieal 
meehanisms. We compare SFA results with the full numerical solutions of the time-dependent 
Schrodinger equation (TDSE) within the few-cycle pulse regime. Good agreement between our SFA 
and TDSE model is found for the ATI spectrum. Our model captures also the strong dependence 
of the photoelectron spectra on the carrier envelope phase of the laser field. 
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I. INTRODUCTION 


During the last three decades, advances in laser technology and the understanding of non¬ 
linear processes in laser-matter interactions have led to production of few-cycle femtosecond 
(1 fs = 10“^® s) laser pulses in the visible and mid-infrared regimes pQE]. By focusing such 
ultrashort laser pulses on a gas target, the atoms are subjected to an ultra-intense electric 
held, with peak held strengths approaching the binding held inside the atoms themselves. 
Such helds are commonly used as a tool to explore the interaction between strong electromag¬ 
netic coherent radiation and an atomic or molecular system with unprecedented spatial and 
temporal resolution [3]. Phenomena such as high-order harmonic generation (HHG) (HE], 
above-threshold ionization multi-photon ionization and multi-electron effects laiH] , are 

routinely studied. These effects can be used to generate attosecond pulses in the extreme 
ultraviolet PIlo] or even soft X-ray regime m- They can also be used to extract either 
information about the laser pulse electric held itself [I2] , or about the structure of the target 
atom or molecule [I31II3]. 

Since electronic motion is governed by the waveform of the laser electric held, an impor¬ 
tant quantity to describe the electric held shape is the so-called absolute phase or carrier-to- 
envelope-phase (CEP). Control over the CEP is paramount for extracting information about 
electron dynamics, and to retrieve structural information from atoms and molecules dSHn]. 
Eor instance, in HHG an electron is liberated from an atom or molecule through ionization, 
which occurs close to the maximum of the electric held. Within the oscillating held, the elec¬ 
tron can thus accelerate along oscillating trajectories, which may result in re-collision with 
the parent ion, roughly when the laser held approaches a zero value. Control over the CEP 
is particularly important for HHG, when targets are driven by laser pulses comprising only 
one or two optical cycles. In such situation CEP determines the relevant electron trajecto¬ 
ries, i.e. the CEP determines whether emission results in a single or in multiple attosecond 
bursts of radiation [ElIIl]. 

The inhuence of the CEP on electron emission was also demonstrated in an anti¬ 
correlation experiment, in which the number of ATI electrons emitted in opposite directions 
was measured da [IS]. Since the first proof of principle experiment im. the stereo ATI 
technique has established itself as a direct measure of the CEP, and demonstrated its 
ability for single shot measurements even at multiple kHz laser repetition rates. The sen- 
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sitivity to the CEP arises from contributions of both, bound-free and the re-scattering 
continuum-continuum transitions of the atomic or molecular target, which are embedded in 
the photoelectron distribution of ATI [Hj. Hence, this mechanism can be used to extract 
structural information about the target atom or molecule. 

Laser induced electron diffraction (LIED) was suggested early on as a technique that 
uses the doubly differential elastic scattering cross section to extract structural information 
[20H22]- Meeting the requirements to extract structural information has, however, proven 
difficult due to the stringent prerequisites on the laser parameters. During recent years, 
the development of new laser sources has dramatically advanced, leading to first demon¬ 
strations [I11I23H26], and the successful retrieval of the bond distances in simple diatomic 
molecules with fixed-angle broadband electron scattering j22]- Recently, Pullen et al. m 
have exploited the full double differential cross section to image the entire structure of a 
polyatomic molecule for the first time. An important next step to exploit the full potential of 
the re-collision physics is the exploitation of the intrinsic time resolution of LIED to extract 
dynamic structural information. The key for such a goal is, however, a comprehensive and 
complete understanding of the ATI process and its theoretical description [191 127II32] . 

The aim of our paper is to revisit the strong held approximation model of M. Lewenstein 
for ATI for few-cycle infrared (IR) laser pulses and to compare it with the numerical solution 
of the TDSE in one (ID) and two (2D) spatial dimensions for an atomic system [29]. For 
simplicity, our analytical atomic model is based on a non-local potential, which can be 
considered a short-range (SR) potential. In order to verify the validity of our analytical 
SR model, and to understand how its predictions compare with a true Coulomb potential, 
we numerically integrate the TDSE for the hydrogen atom and compute the photoelectron 
energy and momentum distributions. 

This article is organized as follows. In Section II, we review the theory which describes the 
ATI process within the Strong Field Approximation; in particular, we present the derivation 
of the transition amplitude for both the direct and re-scattered electrons. We develop in 
detail the mathematical foundations towards the final results by starting from the Hamilto¬ 
nian, which describes the atomic system and the TDSE associated to it. In Section HI, we 
introduce the model for our atomic system, that uses a particular form of a non-local short- 
range potential. The matrix elements to describe the ionization and re-scattering processes 
are then computed. In section IV, the ATI spectra for the ID line and 2D case are numeri- 
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cally calculated and compared with numerical results obtained from TDSE calculations. In 
addition, we discuss the effect of the CEP on the spectra, calculated from our analytical SEA 
model. Einally, in section V, we summarize the main ideas and present our conclusions. We 
give an outlook on extending this analytical model to more complex atomic and molecular 
systems. 


II. STRONG FIELD APPROXIMATION: TRANSITION PROBABILITY AMPLI¬ 
TUDES 

The interaction of a strong electric field with an atomic or molecular system is described 
within Quantum Mechanics by the time dependent Schrodinger equation that captures both 
the evolution of the (electronic) wave function and the time evolution of the physical observ¬ 
ables. The numerical solution of the TDSE offers a full quantum mechanical description of 
the laser-matter interaction processes, and has been extensively used to study several phe¬ 
nomena, such as HHG [3311^ and ATI |36l - Hn] in atomic and molecular systems. However, 
the full numerical integration of the TDSE in all the degrees of freedom of the system is often 
a laborious, and sometimes an impossible task to perform from a numerical and computa¬ 
tional points of view. Moreover, a physical interpretation of the numerical TDSE results 
and the extraction of information from the time evolved wave function is highly nontrivial 
for an ab initio technique. 

Hence, from a purely theoretical point of view, it would be desirable to solve the TDSE 
analytically for the ionization process. This is one of the main steps in all laser-matter 
interaction phenomena, and it represents a formidable and challenging assignment. Here, 
we discuss an alternative method to calculate photoelectron spectra from atomic systems 
by analytically solving the TDSE under the so-called Strong Field Approximation. This ap¬ 
proach dates back to Keldysh [H] , and has since been employed by many other authors [271- 
[221112H1S] • It is worth noted that SEA provides a quantum framework and extension of the, 
so called, “simple man’s” or “three step” or “re-collision” model, usually attributed to P. 
Corkum ra. K. Kulander [481 [49] H. Muller (cf. [50] for an extensive review; for earlier 
quantum formulation of “Atomic Antennas” see Ref. 0; for other pioneering contributions 
see Refs. [S2HS1]). 
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Ionization driven by strong fields 


Let us consider an atom under the influence of an ultra-intense laser field. In the limit 
when the wavelength of the laser, Aq, is larger compared with the Bohr radius, Uq (5.29x 10“^^ 
m), the electric field of the laser beam around the interaction region can be considered spa¬ 
tially homogeneous. Consequentially, the interacting atoms will not experience the spatial 
dependence of the laser electric held and, hence, only its time-variation is taken into account. 
This is the so-called dipole approximation. In this approximation, the laser electric held can 
be written as: 

E{t) = So f{t)sm{uJot + (f)o)ez. (1) 


The field of Eq. ( 1 ) has a carrier frequency uq = where c is the speed of light. So the 
field peak amplitude or strength, and we consider that the laser held is linearly polarized 
along the z-direction. f{t) denotes the envelope of the laser pulse and the parameter 
defines the CEP. 

The TDSE is defined (atomic units are used throughout this paper unless otherwise 
stated) by: 

d 

( 2 ) 




where the Hamiltonian operator, H^ describes the laser-atom system and is the sum of two 
terms, i.e. 


H = Ho + U, (3) 

where, //q? is the so-called laser-free Hamiltonian of the atomic or molecular system 

H« = ~ + V(r), (4) 

with E(r) the atomic or molecular potential, and U = — gE(t) • r, is the dipole coupling 
which describes the interaction of the atomic or molecular system with the laser radiation, 
written in the length gauge and under the dipole approximation. Note that in atomic units, 
the electron charge, denoted by g, is g = — 1 a.u. 

We shall restrict our model to the low ionization regime, where the SEA is valid [271- 
EH HU sal is] and successfully describes the laser-matter interaction processes. Therefore, 
we consider the strong field or tunneling regime, where the Keldysh parameter 7 ' = y^Ip/2Up 
[Ip denotes the ionization potential of the atomic or molecular system and, Up = the 

etu/o 
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ponderomotive energy taken by the electron in the oscillating electromagnetic field) is less 
than one, i.e. 7 ' < 1 . In addition, we assume that the remaining Coulomb potential, C(r), 
does not play an important role in the electron dynamics once the electron appears in the 
continuum. These observations, and the following three statements, define the standard 
SFA, namely: 

(i) The strong field laser does not couple with any other bound state. This means that 
only the ground state, | 0 ), and the continuum states, |v), are taken into account in the 
interaction process; 

(ii) There is no depletion of the ground state, i.e. the ponderomotive energy is lower than 
the saturation energy of the system {Up < Ugat)', Despite this assumption, including 
depletion effects, e.g. by including ionization rates according to the Ammosov-Delone- 
Krainov theory (ADK rates [l 2 ]) provides no a particular challenge; 

(in) The continuum states are approximated by Volkov states; more precisely (cf. [221 SSI 
SS]) the continuum-continuum matrix elements are decomposed in the basis of scat¬ 
tering states, corresponding to waves with a fixed outgoing (kinetic) momentum Pe, 
into the most singular part and the rest, which is treated then in a perturbative man¬ 
ner [29]. In such decomposition, the most singular part corresponds exactly to approx¬ 
imating the scattering states by plane waves, i.e. Volkov solutions. Corrections with 
respect to the less singular part of the continuum-continuum matrix elements describe 
re-scattering and re-collision events. 

Based on the statement (i), we propose a state, |5'(t)), that describes the time-evolution 
of the system by a coherent superposition of the ground, | 0 ), and the continuum states, 

|v) [221112]: 

|4'(i)>=e"»‘(a(t)|0>+ y'rf=v(.(v.t)|v)). (6) 

The factor, a{t), represents the amplitude of the ground state which will be considered 
constant in time, a{t) ~ 1 , under the assumption that there is no depletion of the ground 
state. The last step follows directly from statement (ii). The pre-factor, represents 
the phase oscillations which describes the accumulated electron energy in the ground state 
{Ip = —Eq is the ionization potential of the atomic system, with Eq, the ground state energy 
of the atomic system). Furthermore, the transition amplitude to the continuum states is 
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denoted by 6 (v, t) and it depends both on the kinetic momentum of the outgoing electron 
and the laser pulse. Therefore, our main task will be to derive a general expression for the 
amplitude b{v,t). In order to do so, we substitute Eq. ([^ in Eq. ([^ and by considering, 
Ho\0) = —/p|0), and [—+ E(r)] |v) = ^|v), the evolution of the transition amplitude 
becomes: 

i j d\b{y,t)\w) = J d\ (^ + Ip^b{v,t)\v) + E{t) ■ r\0) 

+E(t) • J d^v [zVv&(v, t) + 6 (v, t)r] |v). ( 6 ) 

Note that we have assumed that the electron-nucleus interaction is neglected once the 
electron appears at the continuum, i.e. E(r)|v) = 0 , which corresponds to the statement 
(iii). Therefore, by multiplying Eq. ([^ by (v'| and after some algebra, the time variation of 
the transition amplitude 6 (v, t) reads: 

^(v, t) = -i + Ip^ 6 (v, t) + iE{t) ■ d(v) 

E(t) • Vv 6 (v, t) - iE{t) ■ J dV b{V, t)g(v, v'). (7) 

The first term on the right-hand of Eq. Q represents the phase evolution of the electron 
in the oscillating laser field. In the second term we have defined the bound-free transition 
dipole matrix element as: 

-{v|r|0> = d(v), (8) 

and finally, the last two terms describe the continuum-continuum transition, Vvb{v,t), 
without the influence of the scattering center, and by considering the core potential, 
J d^v' 6 (v', t)g(v, v'). Here, g(v,v'), denotes the re-scattering transition matrix element, 
where the potential core plays an essential role: 

(v|r|v') = g(v,v'). (9) 

In the following, we shall describe how it is possible to compute the transition amplitude, 
6 (v,t), by applying the zeroth and first order perturbation theory to the solution of the 
partial differential equation Eq. 0 . Therefore, according to this perturbation theory, we 
split the solution of the transition amplitude, b{v,t), into two parts: 6 o(v,t) and 6 i(v,t), 
i.e. b{v,t) = 6 o(v,t) -|- 6 i(v,t). The zeroth order of our perturbation theory 6 o(v,t) will be 
called the direct term. It describes the transition amplitude for a laser-ionized electron that 
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will never re-scatter with the remaining ion-core. On the other hand, the first order term, 
named re-scattered term, bi(v,t), is referred to the electron that, once ionized, will have a 
certain probability of re-scattering with the potential ion-core. 


Direct transition amplitude 


Let us consider the process where the electron is ionized without probability to return 
to its parent ion. This process is modeled by the direct photoelectron transition amplitude 
bo{v,t). As the direct ionization process should have a larger probability compared with 
the re-scattering one [20], one can neglect the last term in Eq. Q. This is what we refer to 
zeroth order solution: 


dtbQ{v,t) 



6o(v,t) -F iE{t) ■ d(v) -F E{t) ■ Vv6o(v,t). 


( 10 ) 


The above equation is a first-order inhomogeneous differential equation, which is easily 
solved by conventional integration methods (see e.g. [55]). Therefore, the solution can be 
written as: 


f,o(v, t) =i [ di! E(t') • d (v - A{t) + A(t')) 
Jo 


X exp ( J' df[{v - A{t) + A(t"))72 + 4] ) • 


( 11 ) 


Here, we have considered that the electron appears in the continuum with kinetic mo¬ 
mentum v(t') = V — A{t) -|- A(t') at the time t', where v is the final kinetic momentum (note 
that in virtue of using atomic units, where the electron mass m = 1, the kinetic electron 
momentum Pe and the electron velocity v have the same magnitude and direction), and 
A(t) = — /* E{t')dt' is the vector potential of the electromagnetic field. In particular, the 
vector potential at the time when the electron appears at the continuum t' is denoted by 
A(t') and at a certain detection time t, the vector potential reads A(t). In addition, it is pos¬ 
sible to write Eq. 0 as a function of the canonical momentum p, defined by p = v — A(t), 
and therefore the probability transition amplitude for the direct electrons simplifies to m-- 


{p,t) =i [ dl!E{t') ■d{p + A{t')) 
Jo 


X exp (-^ / di[{p +A{i))‘^/2 +Ip]\ . 


( 12 ) 


This expression is understood as the sum of all the ionization events which occur from 
the time f to t. Then, the instantaneous transition probability amplitude of an electron 
at a time f, at which it appears into the continuum with momentum v(t') = p + A(t'), is 
dehned by the argument of the integral in Eq. (@1. Furthermore, the exponent phase factor 
in Eq. (12) denotes the “semi-classical action”, S'(p,t,t'), that defines a possible electron 
trajectory from the birth time t' until the “detection” one t [29]: 

rt 


S(p,t,i')=/ df [(p +A(f)f/2 + 4] 


(13) 


As our purpose is to obtain the hnal transition amplitude 6o(p, t), the time t will be hxed 
at the end of the laser field, t = tp. For our calculations, we shall define the integration time 
window as: t: [0,tF]- Therefore, we set, E(0) — E{tp) = 0, in such a way to make sure that 
the electromagnetic field is a time oscillating wave and does not have static components. 
The same arguments are applied to the vector potential A(t). We have defined the laser 
pulse envelope as f{t) = sin^(^) where W denotes the number of total cycles. 


Re-scattering transition amplitude 


In order to find a solution for the transition amplitude of the re-scattered photoelectrons, 
6i(v,t), we have considered the re-scattering core matrix element g(v, v') term of Eq. 0 
different than zero, i.e. g(v, v') ^ 0. In addition, the first-order perturbation theory is 
applied to obtain 6i(v,t) by inserting the zeroth-order solution 6o(p,t) in the right-hand 
side of Eq. Q. Then, we obtain bi{p,t) as a function of the canonical momentum p as 
follows: 

-i,S(p,t,t')]E(t')- [ df 
Jo 

X E(t") • d (p' + A(t")) exp [-iSip’, t', t'% 


/dVg(p + A(t'),p' + A(0) ^ ^ 

J (14) 


bi{p,t) = — / dt'exp\ 


This last equation contains all the information about the re-scattering process. In partic¬ 
ular, it is referred to the probability amplitude of an emitted electron at the time t", with an 
amplitude given by E(t") ■ d (p' -|- A(t")). In this step the electron has a kinetic momentum 
of v'(t") = p' + A(t"). The last factor, exp [— z5'(p', t', t")], is the accumulated phase of an 
electron born at the time t" until it re-scatters at time t'. The term, g(p -|- A(t'), p' -|- A(t')), 
contains the structural matrix element of the transition continuum-continuum at the re- 
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scattering time t'. At this particular moment in time, the electron changes its momentum 
from p' + A(t') to p + A(t'). We stress out, however, that the term g(v, v') does not neces¬ 
sarily imply that the electron returns to the ion core. In addition, exp [—i5'(p, t, t')] defines 
the accumulated phase of the electron after the re-scattering from the time t' to the “final” 
one t when the electron is “measured” at the detector with momentum p. In particular, 
note that the photoelectron spectra, |6(p, is a coherent superposition of both solutions, 
bo(j),tF) and bi{p,tF): 

|&(p,^f)P = | 6 o(p,^f) + &i(p,tF)P, 

= |6o(p,tF)P + \bi{p,tF)\‘^ + bo{p,tF)bl{p,tF) + c.c. (15) 


So far we have formulated a model, which describes the photoionization process leading 
to two main terms, namely, a direct 6 o(P)^f) and a re-scattering 6i(p,tF) one. As the 


complex transition amplitude, Eq. (12), is a “single time integral”, it can be integrated 
numerically without major problems. However, the multiple time (“2D”) and momentum 


(“3D”) integrals of the re-scattering term, Eq. (14), present a very difficult and demanding 
task from a computational perspective. In order to reduce the computational difficulties, 
and to obtain a physical meaning of the ATI process, we shall employ the stationary phase 
method to evaluate these highly oscillatory integrals. 

The fast oscillations of the momentum p' integral for the electron re-scattering transition 
amplitude, bi{p,t), suggests to use the stationary-phase approximation or the saddle point 


method to solve Eq. (14). This method is expected to be accurate, when both the Up and 
the Ip, as well as the involved momentum v and v', are large. As the quasi-classical action 
S{p',t 'is proportional to Ip, Up and v'^, the phase factor, eyip{—iS{p',t',t"), oscillates 


very rapidly. Then, the integral over the momentum p' of Eq. (14) tends towards zero except 
near the extremal points of the phase, i.e. Vp/5'(p') = 0. Thus, the main contributions to 
the momentum integral are dominated by momenta, p(, which satisfy the solution of the 
equation: Vp/<S'(p')|p/^ = 0. These saddle point momenta read: 


P. 


1 

r 


h" 


A(i)di. 


(16) 


Here, r = t' — t" is the excursion time of the electron in the continuum. In terms 
of Classical Mechanics, these momenta roots p( are those corresponding to the classical 
electron trajectories because the momentum gradient of the action can be understood as 
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the displacement of a particle [SE]- As the momentum gradient of the action is null Ar = 
Vp/5'(p', t', t") = 0, the considered electron trajectories, r(t), are for an electron that is born 
at the time t" at a certain position r(t") = Tq. Then, after some time t' the electron returns 
to the initial position r(t') = Tq with an average momentum p'. 

Therefore, the function 5'(p', t', t") can be expanded in Taylor series around the roots p^ 
and the transition amplitude for the re-scattering electrons bi{p,t) becomes: 


6 ,(p, t)=- f *[(p+A(t))V2+4] E(t') . f dfg (p + A(t'), p' + A(t')) 


X 


TT 




E(i") ■ d (p^ + A(t")) e-'-f*" [(p;+A(i))’/2+/p| 


(17) 


Here, we have introduced a smoothing parameter, e, to avoid the divergence at the time 


t' = t". Note that the 3D momentum integral on p' of Eq. ( [14[ ) can then be solved by: 

j o^p7(pl)exp 


A(p;) + ^Vj,A(p') 


(P' - Ps? 


7T 


s + 




/(P'J- (18) 


With the last equation we have substantially reduced the dimensionality of the problem 
from a 5D integral to a 2D integral. As the computing time depends on the dimensionality 
of the integration problem, this reduction is extremely advantageous from a computational 
viewpoint. Moreover, with the saddle point method a quasi-classical picture for the re¬ 
scattering transition amplitude is obtained similar to the approach described in 1221 EZ]. 

The main problem to calculate the ATI spectrum is then the computation of the bound- 
free transition dipole matrix element, d(v), and the continuum-continuum transition re¬ 
scattering matrix element g(v, v') for a given atomic system. In the next section, we shall 
introduce a short-range potential model in order to compute the transition matrix elements 
and the final photoelectron momentum distribution analytically. 


III. ABOVE-THRESHOLD IONIZATION IN ATOMIC SYSTEMS 

In this section, as a test case for our model, we chose a non-local atomic potential with 
the purpose of computing both the direct and the re-scattering transition amplitudes. These 
terms involve the dipole and the continuum-continuum matrix elements defined by Eqs. (|^ 
and Then, our main task will be devoted to analytically find the wavefunctions for the 
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ground and scattering states of our test potential. The Hamiltonian, H{p, p'), of the atomic 
system in the momentum representation can be written as: 

-^(p,p') = Y^(p-p')+ ''^(P>P')> (19) 


where the first term on the right-hand side is the kinetic energy operator, and the second 
one is the interacting non-local potential ^(p, p^). By using such Hamiltonian, we write the 
stationary Schrodinger equation as follows: 

H{p,p')^{p) = I (fp'H{p,p')^{p'), 

^^(p) = Y y'^ p''^(P“P0^(p')- 70(p) y cf^pV(p')^(p')> (20) 

where E denotes the energy of the wavefunction T(p). Note that we have defined the 
non-local potential as r(p,p') = — 70(p) (f>{p'), which describes the attraction between the 
electron and the nucleus [22] • This potential has been chosen such that it assures analytical 
solutions of the continuum or scattering states, i.e. for states with energies E > 0. Note 
that the ground state can also be calculated analytically. 7 is understood as a screening 
parameter and 0(p) is an auxiliary function defined by: 


'#'(p) = 


1 


( 21 ) 


where the parameter T is a constant related with the shape of the ground state. In order to 
analytically obtain the ground state, d^o(p)) we solve the stationary Schrodinger equation 
in the momentum representation: 

-2 (fp'-qtoip') 


— dfo(p)- - — 

2 ^/^+f^J yp'2^p2 


= £^o^'o(p), 


( 22 ) 


where the parameter 7 is related to the ionization potential, E, of the atomic species under 


study. To solve Eq. (22) we consider (p = f 


-v/p' 2 +r 2 


as a new parameter and write the 


eigenenergy Eq = —Ip. Therefore, the final solution reads: 


fo)?) = 


v 


(23) 


\/(?Tr7(? + Ip) 

where, J\f = ^ (p denotes a normalization constant. Dividing the last formula by -|- T^ 
and taking the volume integral on p, we obtain: 

(^p 




(p2 + r2)(2f+ 4)' 


(24) 
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The solution of the last integral in Eq. (24) gives us the relation between the parameters 


Ip, r and 7 : 


7 / dU d(psm(r> -^- 

Jo Jo Jo (p2 + r2)(f + 4) 


7 


r + ■\/2ip 


= 1, 


= 1. 


(25) 


This formula allows us to control the parameters T or 7 , in such a way as to match the 
Ip of the atomic system. Furthermore, by using the normalization condition for the bound 
states, we calculate the normalization constant, Af, as well as the analytical ground wave 
function To(p). This normalization factor reads: 





(26) 


So far we have obtained, analytically, the ground state of our non-local potential model. 
This ground state will allow us to calculate the bound-free transition dipole matrix element 
by using Eq. ([^. The free or continuum state is approximated as a plane wave of a given 
momentum, po, and therefore the bound-free transition dipole matrix in the momentum 
representation reads: 

d(Po) = * J ^o(p')Vp'(5(p' - Po) dV- (27) 

By employing properties of the Dirac delta distribution, d(po) is computed via d(po) = 
— ^Vp/To(p') Ipjj- After some elementary algebra, we obtain the transition dipole matrix : 


d(Po) 




_ K _^ 

iP'^+nk^+ip) J 


= iAfpo 


(Po+r^)+(^+/p) 

(p§+r2)2(f+/p)2' 


5 

PO 


(28) 


The second important quantity to be calculated before evaluating the whole transition 
amplitude 6(p,t) is the transition continuum-continuum matrix element, g(p,pO. Hence, 
we need to find the scattering or continuum wave functions of our model potential. Next, we 
shall calculate the scattering states by analytically solving the time independent Schrodinger 
equation in the momentum representation for positive energies. 
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Scattering waves and continuum-continuum transition matrix element 


Let us consider the scattering wave, (p), with asymptotic momentum po, as a coherent 
superposition of a plane wave and an extra correction (5d'pg(p): 

^Po(p) = '^(p-Po) + <^^Po(p)- (29) 

This state has an energy E = Po/2. Then, the Schrodinger equation in momentum 
representation reads: 


d^p' d'po(p) 


f ^P.(P) = I ^p.(P) - j 


— - ^ 1 hT 


(P) = 


7 


+ 


7 


C^P' <^^po(p) 

+ lVpo + r" ' VPTpy x/p'2 + r2 ' 


/ 


2 2 ' ” ~^ / n . -r-.n / O . -r-io ' / O . / / .o ■ (^^) 

To analytically solve the last equation, we apply elementary algebra and the following 
Dirac delta distribution properties: ^ ^ _ n „—i r p '5(p-po) _ i 

Finally, the correction ^Tp^, results: 

B{po) 


) 5(p — Po) = 0 , and f _ 


<^^Po(P) = 


^e 


(31) 


Vp^ + r^jpg-p* 

Here, e, is a smooth parameter to avoid the divergence at p = po and -B(po) is a constant, 
which depends on the asymptotic momentum po- The constant -B(po) is defined by: 


5(Po) = -27 


1 


LvS+T^ 




(32) 


where (f' = f order to obtain i?(po), we proceed analogously to our procedure 


for obtaining Eq. (|24|). Consequently, for Eq. 

:toK 
27 


the left and right-hand sides. These factors cancel each other and the constant H(po) reads: 

g(Po)= , , ’■ (33) 


we obtain the same quantity, (p', on both 
3 h othe 
471^^7 


(Po + r2)2V bol + 

Einally, the scattering wave functions can be written as: 

i?(Po) 


^Po(P) = <^(P - Po) 


\/^ 


(34) 


+ fe 

The later equation tells us that the correction, h4'p(,(p), to the plane wave is a function 
of the parameters of the atomic potential, T and 7 . Therefore, the re-scattering process 
will depend on the shape of the potential. However, in the limit when the momentum po 
goes to infinity this correction term vanishes, i.e. limpQ_j.oo <^fl^po(p) = 0, and then the atomic 
potential does not play any role in the re-scattering process. 
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Continuum-continuum transition matrix element 


Let us consider the scattering waves obtained in Eq. (34) and evaluate the continuum- 
continuum transition matrix element of Eq. I§. i.e. 

/ + (X) 

^pl(p) Vp d'p2(p)c^P, 

-oo 




—oo 
+ 00 


X Vr 


<^(p - Pi) 
<^(P-P 2 ) + 




-\- r^) 2 {p\ — 
5(P2) 


d^p. 


(35) 


ip^ + r^) 2 (p2 — p2 _|_ 

As the first-order perturbation theory has been considered along our derivations, all 
quadratic or superior terms in 7 , e.g. i?*(p^)i?(p 2 ), are neglected. Therefore, we obtain: 

1 


/ H-oo 

<^(p -Pi)V 

-00 


+*i^*(Pi) 


p+00 


-i- r^) 2 (p2 — p2 _ 

1 


d^p 


Vpd(p-P2)d p. 


(36) 


'-00 (p^ -h r^) 2 (p2 _ p2 _|_ ^g^ 

The last momentum integrals are solved by applying the same Dirac delta distribution prop¬ 


erty used in Eq. (28) and, after some algebra, the transition matrix continuum-continuum 


element for our model potential reads: 


g(Pl>P2) = *^(P2)Pl 


?>p\-p^+2T+ie 


- *-B*(Pi)P 2 


3p|—p^+ 2 r— 


(37) 


_ (pf +r2 ) 2 (p2 _p 2 _ 2 J [ (p| -|-r2 ) 2 (p2 _p 2 _|_ 2 ^ 

At this point, we have obtained all the required elements to evaluate both the direct and 
the re-scattering transition amplitude terms defined according to Eqs. (12) and ( [1^ . The 
developed model is an alternative way to describe the ATI process mediated by a strong 
laser pulse. The method is physically intuitive, and can be understood on the basis of a 
quasi-classical picture, i.e. electron trajectories. This is the main difference of our approach 
in comparison to the numerical solution of the TDSE, whose physical interpretation is, in 
spite of its accuracy, frequently challenging. The main advantage of the proposed model is 
that Eqs. (12) and 10 give a clear physical understanding of the ATI process and provide 
rich information about both the laser field and the atomic target which are encoded into 
the complex transition amplitude 5(p,t) = 6o(P;^) + 5i(p,t). The exact analytical solutions 
of those direct and re-scattering transition amplitudes are however not nontrivial to obtain 
if no approximations are considered. In particular, for the re-scattering photoelectrons, the 
solution is even more complex and depends, generally, of the laser electric field shape. 
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In the next Section, we numerically integrate both terms, i.e. 6o(p, t) and 6i(p, t), for the 
non-local potential and compare those results to the numerical solution of the TDSE. 


IV. RESULTS AND DISCUSSION 


The numerical integration of Eqs. ( 12 ) and 0 has been performed by employing a 
rectangular rule with dedicated emphasis on the convergence of the results. As the final 


momentum distribution Eq. (15) is “locally” independent of the momentum p, i.e. |6(p,t)p 
can be computed concurrently for a given set of p values, we have optimized the calculation 
of the whole transition amplitude, |6(p, t)p, by using the OpenMP parallel package [HH]- The 
final momentum photoelectron distribution, |6(p, t)p, is computed both in a ID-momentum 
line along p^, and in a 2 D-momentum plane {py, Pz)- We shall compare these results with the 
numerical solutions of the TDSE in one (ID) and two (2D) spatial dimensions, respectively. 

In case of the ID calculations for the ATI spectra, the momentum grid was symmetrically 
defined with a length of = 4.0 a.u., and a step size of 5pz = 0.02 a.u. The parameters 
of the non-local potential are fixed to E = 1 and 7 = 38 a.u., in such a way as to match 
the ionization potential of the hydrogen atom. Ip = 0.5 a.u. Note that several values of T 
and 7 can be employed to obtain the same Ip. Therefore, these parameters are chosen to 


match the ground state wave function, Eq. (23), of our SR potential model with the shape 


of the ground state wave function of an actual hydrogen atom. We use in our simulations 
an ultrashort laser pulse with central frequency ujq = 0.057 a.u. (wavelength A = 800 nm, 
photon energy, 1.55 eV), peak intensity /q = 1 x 10^^ W -cm"^, with a sin^ envelope shape 
with Vc = 4 total cycles (this corresponds to a full-width at half-maximum FWHM = 2.67 
fs) and a CEP 0o = 0 rad. The time step is fixed to 5t = 0.2 a.u., and the numerical 
integration time window is t\ [Ojtp], where tp = 11 fs and Tq = 2'n/ujQ denote the 

final “detection” time and the cycle period of the laser field, respectively. 

Figure [T] shows the final photoelectron distribution or ATI spectra, in logarithmic scale, 
as a function of the ponderomotive energy, Up, for electrons with positive momenta along 
the p^-direction. Fig. [^a) depicts the total contribution, Eq. (15), meanwhile Fig. ib) 
shows the contribution of both the direct |6o(p,t)P and re-scattering terms |6i(p,t)p. For 
completeness, the interference term, 6 int(pU) = ^o(PT) ^i(pU) + ^o(pU) ^i(pT) E included 
as an inset of Fig. [^a). The first clear observation is that each term contributes to different 


16 






regions of the photoelectron spectra, i.e. for electron energies Ep^ < 3Up the direct term 
l^o(P) t)\^ dominates the spectrum and, on the contrary, it is the re-scattering term, |6i(p, 
the one that prevails in the high-energy electron region. In addition, we observe that the 
interference term follows the trend of the direct one (see the inset of Fig. [^a)) and does 
not play any role for electron energies Ep^ > bUp. We shall see next that both direct and 
re-scattering terms are needed in order to adequately describe the ATI process. 



Energy (Up) 


(a) Full ATI photoelectron spectrum and 
interference term (inset) 



(b) Direct and re-scattering photoelectron 
spectra 


FIG. 1: (color online) Photoelectron ATI spectra (in logarithmic scale) as a function of the 
ponderomotive electron energy Up computed by using our quasi-classical model and for 


each one of the transition terms: (a) Total photoelectron spectra, Eq. (15), (red line) with 
the interference term in the inset (magenta line), (b) Direct photoelectron spectrum 
|&o(Pz)^f)P and re-scattering photoelectron spectrum \bi{pzU'F)\^ are depicted in blue solid 
and black dashed lines, respectively. The vertical dashed lines correspond to the classical 

2Up and IQUp cutoffs (see the text for details). 


To confirm that our model is able to capture the left-right asymmetry [TH], in Fig. 
we compute ATI spectra for electrons with positive and negative momenta along the p^- 
direction. Fig.j^a) shows the results of our quasi-classical model, meanwhile in Fig. [^b) the 
TDSE in ID is used. The photoelectrons with negative {positive) momentum are convention¬ 
ally named left (right) electrons and correspondingly the photoelectron spectra associated 
are labeled by |5 l(Pz)0o)P and \b^^(pz,(l)o)\^, respectively. 

The photoelectron spectra computed by using the numerical solution of the TDSE in ID, 
Eig. [^b), allow us to evaluate the accuracy of our quasi-classical ATI model. The numerical 
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(a) SFA Calculations (b) TDSE Calculations 

FIG. 2: (color online) Comparison of the ATI spectra for an hydrogen atom, (a) 
Photoelectron energy distribution (in logarithmic scale) for the emitted electrons with 
negative (green dark line) and positive (red line) momentum obtained by the integration of 
our derived full transition amplitude \h{pz^tY)\^. (b) The same as in (a) but computed by 
the numerical solution of the TDSE in ID. The vertical dashed lines correspond to the 
classical 2Up and lODp cutoffs (see the text for details). 


integration of the TDSE is performed by using the Split-Spectral Operator algorithm [59] 

and we use the EFTW iOI to evaluate the kinetic energy operator of our Hamiltonian 

^ -^2 ^ 

H = ^ + V{z) + zE(t). For the present numerical solution of the TDSE, we have fixed the 
position grid step to 5z = 0.2 a.u., with a total number of points = 17000. The ground 
state is computed via imaginary time propagation with a time step of 5t = — 0.02 i and the 
soft-core Coulomb potential is given by: V{z) = — ^/A_ . The parameter a — 2 a.u., is 
chosen in such that the ground state yields the ionization potential of the hydrogen atom, 
i.e. Ip = 0.5 a.u. 

The strong-field laser-matter interaction is simulated by evolving the ground state wave 
function in real time, with a time step of 5t = 0.02 a.u., and under the action of both the 
atomic potential and the laser field. The laser pulse parameters are the same as those used 
to compute the results of Eig. [T} At the end of the laser field tp, when the electric field is 
zero, we compute the final photoelectron energy-momentum distribution |6 tdse(PzTf)P; by 
projecting the “free” electron wave packet, Tc(.2,tF), over plane waves. The wave packet 
d^c(^) ^f), is computed by smoothly masking the bound states from the entire wave function 
function d^(.2,tF) via: Tc(2:,tF) = h(z)^{z,tp), where, h(z) = exp(—(^^)^) is a gaussian 
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filter. 


Figure demonstrates good qualitative agreement between the photoelectron spectra 
calculated with our quasi-classical model and those obtained by the numerical solution of 
the TDSE in ID. The left-right photoelectron spectra show the expected two cutoffs defined 
by 2Up and IQUp (black dashed fines) which are present in the ATI process [HI EH]- This 
shows that our approach is a reliable alternative for the calculation of ATI spectra. Our 
model furthermore captures the left-right dependence of the emitted photoelectrons as shown 
in Fig. [^a), and in comparison with the TDSE shown in Eig. |^b). The ability to capture 
this dependence and its features is especially important for applications to methods such as 
LIED which relies on large momentum transfers and backscattered electron distributions. 
For instance, photoelectrons ejected towards the left differ substantially from those emitted 
to the right for the case when a few-cycle driving pulse is used. According to the quasi- 
classical analysis of Section |TI| one can then infers that electron trajectories emitted towards 
the right have larger probability to perform backward re-scattering with the ionic core than 
the electrons emitted towards the left [121 US] . This behavior is clearly reproduced by both 
models shown in Fig. and it is the basis for the stereo ATI technique developed by Paulus 
et al. [T2] . 

Since our model is capable of capturing the general CEP dependence we turn to a more 
detailed investigation on whether our model can reproduce detailed CEP dependence by 
computing the ATI spectra as a function of the absolute laser phase 0o- Henceforth, we 
define the left-right asymmetry A(pz,0o) as visibility: 

|^l(p^,0o)P - |6r(p^,0o)P 


AipzAo) = 


(38) 


\hK{Pz,(t>o)? + |^l(P2,0o)P' 

We compute the ATI spectra for a set of CEP values between 0o = t 180°, and evaluate 
the asymmetry A{pz-,4>o) of Eq. (38). The results are shown in Fig. Our calculated 
asymmetry Aipz, 4 >q) shows a clear dependence on the absolute phase of the laser pulse. 
For instance, when the CEP is 0o = T 90°, the photoelectron spectra show a left-right 
symmetry, which is clearly visible in the energy region between 0 and AUp (see Pig. |^. This 
symmetry can be attributed to the direct term which dominates the photoelectron spectra 
at lower energies and is a consequence of the symmetry of the electric field with respect to 
the envelope maximum. On the other hand, and as we shall see later, the high-energy 
re-scattered electrons do not follow this symmetry. 
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CEP (°) 

FIG. 3: (color online) Asymmetry of the photoelectron energy distribution A{pzi 0o) as a 
function of the CEP. The horizontal dashed white lines denote the 2Up and lOGp cutoffs 
rule for the direct and re-scattering photoelectrons, respectively. The laser pulse and the 
atomic parameters are the same that those used in Fig. [2| 


For the energy range blip < Ep^ < 12Gp, the term |6 lP is less than |6 rP around 0o = 0°. 
This implies that left electron trajectories have less probability to perform backward re¬ 
scattering than those trajectories emitted to the right. Note that this process changes if 
the CEP of the laser pulse is larger than 90° thereby the change is for energies between 
blip < Ep^ < 8Up. In this interval the electron trajectories emitted to the left have a larger 
probability than the ones towards. For low-energy photoelectrons Ep^ < blip, the asymmetry 
oscillates between positive and negative values, which means that the left-right direct pho- 
toelectrons are more difficult to evaluate compared to using re-scattered ones. Thus, these 
results depicted in Fig. [^clearly show that our model describes the typical dependence of 
the ATI spectra on the CEP [121 US] and in particular the backward re-scattering events. 
Our model therefore can be used to describe the absolute phase of the driving IR laser 
pulse. With the purpose to understand the left-right symmetry (or asymmetry) presented 
in Fig. 1^ we compute both the direct and re-scattering terms for two different CEP values, 
namely 0o = 0° and 0o = 90°. The results are depicted in Fig. In the case of 0o = 0° the 
laser pulse is asymmetric with respect to the pulse envelope maximum, i.e. it has a sin(c<;ot) 
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carrier wave. Consequently, from the phase contribution in Eq. (12) of the direct term, one 
can expect that the phase as a function of time is asymmetric as well. It is a consequence 
of the Fourier relation that a temporal asymmetry leads to an asymmetric spectral phase. 
In analogy, the temporal asymmetry of the phase of the direct photoelectron term, leads to 
the final photoelectron momentum distribution being asymmetric with respect 

to the momentum zero axis. This dependence is the origin of the asymmetric shape of the 
left-right emitted photoelectrons shown in Fig. |^a). On the other hand, when the phase of 
Eq. ( [I^ is time symmetric, which is the case of 0o = 90°, i.e. the phase is proportional to 
cos(n;o^)) we infer that the photoelectron spectrum for the direct term should be symmetric. 
This is exactly what we observe in the direct term which is depicted in Fig. |^b). Moreover, 
in both cases the re-scattering term |6i(p^,tF)P is asymmetric with respect to the Pz = 0 
momentum. Hence, from the quasi-classical analysis addressed in Eq. PH and due to the 
few-cycle electric field waveform, the electron trajectories strongly depend on the CEP and 
the left-right momentum asymmetry is visible due to the occurrence and interference of only 
a few emission and re-scattering events. 




(a) 


(b) 


FIG. 4: (color online) Left-right photoelectron momentum distributions (in logarithm 
scale) for the different contributions, direct term |&o(Pz)^f)P (blue dashed line), the 
re-scattering term \bi{pz,t-F)\‘^ (black dashed with points line), and the “total” term, 
\b{pz,tp)\‘^ = \bo{pz,tp) + bi{pz,tF)\'^ (red solid line), to the ATI spectra for two different 
CEP values, 0o = 0°and 0o = 90°are depicted in (a) and (b), respectively. The laser pulse 
and the atomic parameters are the same that those used in Fig. [2} 


In order to complete the analysis, we have extended our numerical calculations of the pho- 
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toelectron momentum distribution for the ATI process from a ID-momentum line to a 2D- 
momentum plane. In Fig. [^we depict results for both models: our analytical quasi-classical 
ATI model (Fig. [^a)) and the exact numerical solution of the TDSE in 2D (Fig. [^b)). We 
find qualitative good agreement between the results of our model and the full numerical 
solution of the TDSE in 2D. We also find that the distribution is symmetric with respect 
to the Py axis. Note that these observations are in good agreement with calculations and 
measurements presented in Refs. [HIHSS]- 



Pz (a.u.) Pz (a.u.) 

(a) SFA calculations (b) TDSE calculations 


EIG. 5: (color online) Comparison between our semi-classical model and the TDSE in 2D 
for an hydrogen system, (a)-(b) Photoelectron ATI spectra |fe(Pz,Py, tF)P (in logarithmic 
scale) computed by employing our model and the TDSE-2D numerical solution, 
respectively. The laser pulse parameters used in these calculations are the same as those 
employed in Fig I Note that the laser field is polarized along the z-direction. 

The comparison shows that our quasi-classical approach can be used to model 2D- 
momentum distributions and even 3D-momentum distributions. However, from the contrast 
between the two models, we infer that our semi-analytical model is limited to photoelectrons 
with high energies. The origin of this discrepancy arises from the approximation made in the 
model with regard to the atomic potential. Statement (hi) relates to the fact that the atomic 
potential is neglected when the electron is born in the continuum. Hence, we expect that 
electrons with lower final energies are not well described by our quasi-classical approach. 

Finally, the main advantage of the analytical model is shown in Fig. which depicts 
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the individual contributions to the 2D ATI spectrum, namely, the direct (Fig. i a)), re¬ 
scattering (Fig. I^b)) and interference term (Fig. |^c)), respectively. The resulting total and 
experimentally accessible ATI momentum spectrum \b{pz,Py,t-F)\^ is shown in Fig. |^d). 



Pz (a.u.) Pz (a.u.) 

(a) Direct Term (b) Rescattering Term 



Pz (a.u.) Pz (a.u.) 

(c) Interference Term (d) Total contribution 


FIG. 6: (color online) Different contributions to the photoelectron spectra for a 
2 D-momentum plane {pz,Py)- ATI photoelectron spectra (in logarithmic scale) as a 
function of the momentum {pz,Py) computed by our quasi-classical model for each term, 
(a) Direct term, (b) re-scattering term, (c) interference term and (d) total contribution. 


The atomic potential and the laser parameters used in these simulations are identical to 
those employed in the calculations for Figs.[T]j^ Analog to the ID calculations, the computed 
photoelectron momentum spectrum for the direct term (Fig. |^a)) shows contributions for 
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electron energies less than 2Up. We find that the contribution of the re-scattering term 
(Fig. I^b)) extends to higher momentum values. Clearly visible is the symmetry of the 
structures about the Py axis for all the terms and a left-right asymmetric shape for electrons 
with P 2 < 0 or P 2 > 0 . 

We like to emphasize the importance of Eqs. and , from which we conclude that 
the form of the calculated ATI spectra depends strongly on the parameters F and 7 of our SR 
potential model. These parameters have a strong influence on the re-scattering term, which 
could get largely suppressed for a particular choice of them. This strong dependence suggests 
that the re-scattering process depends strongly on the atomic target which means that 
particular structural information is encoded in the ATI photoelectron spectra. Consequently 
the proposed semi-analytical model can be used to extract target structure and electron 
dynamics from measured photoelectron spectra. 


V. CONCLUSIONS AND OUTLOOK 

We have studied the photoionization process mediated by a strong laser field interacting 
with an atomic system. We analyzed in detail approximate analytical solutions of the TDSE, 
obtained under the assumption of the strong field approximation, i.e. once an electron is tun¬ 
nel ionized, its dynamics is solely governed by the driving laser, which leads to re-scattering 
or re-collision events. Based on this approach, we have identified and calculated the two 
main contributing terms in the ionization process: the direct and the re-scattering transition 
probability amplitudes. In addition, the bound-free dipole and the re-scattering transition 
matrix elements were analytically computed for a non-local potential. We stressed that this 
is one of the main difference of our developed model that those traditionally found in the lit¬ 
erature for the ATI process. These analytical derivations of the re-scattering matrix element 
allowed us to demonstrate that the re-scattering process strongly depended on the atomic 
target features. A quasi-classical analysis of the re-scattering transition amplitude was per¬ 
formed in terms of the saddle point approximation, which permits linking the dynamics to 
relevant quasi-classical information, i.e. classical electron trajectories. Our analytical results 
suggested that the main contributions to the re-scattering transition amplitude correspond 
to electron trajectories, with significant probability of backward scattering off the ionic core. 

Our model was used to demonstrate that both contributions, the direct and the re- 
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scattering terms, shown left-right asymmetry depending on the carrier envelope phase of 
the laser pulse. This behavior has been confirmed by a comparison with the exact numerical 
solution of the TDSE and we found very qualitative good agreement, particularly in the high 
energy region of the photoelectron spectra. Apart from testing the validity of our model, 
we stress that it presents important advantages, such as the possibility to disentangle the 
effects of both the direct and re-scattered terms. 

We showed also that the model is sensitive to the CEP, and by using the fact that we 
can investigate individual contributions to the photoelectron spectrum, we identified the 
re-scattering term that plays a dominant role by varying its influence based on the atomic 
parameters. These findings confirm that the photoelectron spectra contain structural in¬ 
formation about the re-scattering process, i.e. about the shape of the ground state wave 
function, as well as the bound-free, and free-free matrix transition elements. This depen¬ 
dence implies that atomic structural information can be efficiently extracted with our model 
for methods such as LIED which measure ATI spectra. 

While our aim was the establishment of a basic semi-analytical theoretical framework 
based on the SEA, we note that our approach is applicable to more complex, and thus more 
interesting systems such as molecules. The method could be extended to describe dynam¬ 
ically evolving molecular systems or atomic clusters. It will be interesting to corroborate 
which kind of information could be extracted, and whether one could visualize molecular 
dynamics such as vibrations or dissociation, or if the ground state molecular orbital could 
be reconstructed. We will address these and similar questions in future publications. 


ACKNOWLEDGMENTS 

This work was supported by Ministerio de Economfa y Competitividad through Plan 
Nacional (FIS201T30465-C02-01) and FrOntiers of QUantum Sciences (FOQUS): Atoms, 
Molecules, Photons and Quantum Information (FIS2013-46768-P), the Catalan Agenda de 
Gestio d’Ajuts Universitaris i de Recerca (AGAUR) with SGR 2014-2016, Fundacio Gellex 
Barcelona and funding from LASERLAB-EUROPE, Grant agreement 284464. N.S. was 
supported by the Erasmus Mundus Doctorate Program Europhotonics (Grant No. 159224- 
1-2009-1-FR-ERA MUNDUS-EMJD). A. C. and M. L. acknowledge ERG AdG OSYRIS. 


25 



J.B. acknowledges FIS2014-51478-ERC. 


[1] M. Nisoli, S. De Silvestri, O. Svelto, R. Szipocs, K. Ferencz, C. Spielmann, S. Sartania, and 
F. Krausz, Opt. Lett. 22, 522 (1997). 

[2] M. Hemmer, M. Baudisch, A. Thai, A. Couairon, and J. Biegert, Opt. Exp. 21, 28095 (2013). 

[3] F. Krausz and M. Y. Ivanov, Rev. Mod. Phys. 81, 163 (2009), 

[4] A. L’Huillier and P. Balcou, Phys. Rev. Lett. 70, 774 (1993). 

[5] A. McPherson, G. Gibson, H. Jara, U. Johann, T. S. Luk, 1. A. McIntyre, K. Boyer, and 

C. K. Rhodes, J. Opt. Soc. Am. B 4, 595 (1987). 

[6] P. Agostini, F. Fabre, G. Mainfray, G. Petite, and N. K. Rahman, Phys. Rev. Lett. 42, 1127 
(1979). 

[7] G. Mainfray and C. Manus, Rep. Prog. Phys. 54, 1333 (1991). 

[8] M. Uiberacker, T. Uphues, M. Schultze, A. J. Verhoef, V. Yakovlev, M. F. Kling, J. Rauschen- 
berger, N. M. Kabachnik, H. Schroder, M. Lezius, et a/.. Nature (London) 446, 627 (2007), 

[9] M. Drescher, M. Hentschel, R. Kienberger, G. Tempea, C. Spielmann, G. A. Reider, P. B. 
Corkum, and F. Krausz, Science 291, 1923 (2001). 

[10] K. Zhao, Q. Zhang, M. Chini, Y. Wu, X. Wang, and Z. Chang, Opt. Lett. 37, 3891 (2012). 

[11] F. Silva, S. Teichmann, S. L. Cousin, and J. Biegert, Nature Commun. 6, 6611 (2015). 

[12] G. G. Paulus, F. Grasbon, H. Walther, P. Villoresi, M. Nisoli, S. Stagira, E. Priori, and S. D. 
Silvestri, Nature Phys. 414, 182 (2001). 

[13] C. 1. Blaga, J. Xu, A. D. DiChiara, E. Sistrunk, K. Zhang, P. Agostini, T. A. Miller, L. F. 

DiMauro, and C. D. Lin, Nature (London) 483, 194 (2012). 

[14] M. G. Pullen, B. Wolter, A. T. Le, M. Baudisch, M. Hemmer, A. Senftleben, C. D. Schroter, 
J. Ullrich, R. Moshammer, C. D. Lin, et a/.. Nature Commun. 6, 7262 (2015). 

[15] A. Baltuska, T. Udem, M. Uiberacker, M. Hentschel, E. Goulielmakis, C. Gohle, R. Holzwarth, 
V. S. Yakovlev, A. Scrinzi, T. W. Hansch, et a/.. Nature (London) 421, 611 (2003). 

[16] J. Itatani, J. Levesque, D. Zeidler, H. Niikura, H. Pepin, J. C. Kieffer, P. B. Corkum, and 

D. M. Villeneuve, Nature (London) 432, 867 (2004). 

[17] B. Wolter, C. Lemell, M. Baudisch, M. G. Pullen, X. M. Tong, M. Hemmer, A. Senftleben, 
C. D. Schroter, J. Ullrich, R. Moshammer, et a/., Phys. Rev. A 90, 063424 (2014). 


26 



[18] 1. J. Sola, E. Mevel, L. Elouga, E. Constant, V. Strelkov, L. Poletto, P.Villoresi, E. Benedetti, 
J. P. Caumes, S. Stagira, et al.^ Nature Phys. 2, 319 (2006). 

[19] D. B. Milosevic, G. G. Panins, D. Bauer, and W. Becker, J. Phys. B 39 , R203 (2006). 

[20] T. Zuo, A. D. Bandrauk, and P. B. Corkum, Chem. Phys. Lett. 259, 313 (1996). 

[21] M. Lein, J. P. Marangos, and P. L. Knight, Phys. Rev. A 66, 051404 (2002). 

[22] C. D. Lin, A. T. Le, Z. Chen, T. Morishita, and R. Luccheser, J. Phys. B 43, 122001 (2010). 

[23] J. Xu, C. 1. Blaga, K. Zhang, Y. H. Lai, C. D. Lin, T. A. Miller, P. Agostini, and L. F. 
DiMauro, Nature Commun. 5 , 4635 (2014). 

[24] M. Meckel, D. Comtois, D. Zeidler, A. Staudte, D. Pavicic, H. C. Bandulet, H. Pepin, J. C. 
Kieffer, R. Dorner, D. M. Villeneuve, et al.^ Science 320, 1478 (2008). 

[25] T. Morishita, A. T. Le, Z. Chen, and C. D. Lin, Phys. Rev. Lett. 100, 013903 (2008). 

[26] J. Xu, C. 1. Blaga, A. D. DiChiara, E. Sistrunk, K. Zhang, Z. Chen, A. T. Le, T. Morishita, 
C. D. Lin, P. Agostini, et a/., Phys. Rev. Lett. 109, 233002 (2012). 

[27] F. H. M. Faisal, J. Phys. B 6, 89 (1973). 

[28] H. R. Reiss, Phys. Rev. A 22, 1786 (1980). 

[29] M. Lewenstein, K. C. Kulander, K. J. Schafer, and P. H. Bucksbaum, Phys. Rev. A 51 , 1495 
(1995). 

[30] M. V. Frolov, D. V. Knyazeva, N. L. Manakov, J. W. Geng, L. Y. Peng, and A. F. Starace, 
Phys. Rev. A 89 , 063419 (2014). 

[31] W. Becker, F. Grasbon, R. Kopold, D. B. Milosevic, G. G. Panins, and H. Walther, Adv. At. 
Mol. Opt. Phys. 48, 35 (2002). 

[32] O. 1. Tolstikhin, T. Morishita, and S. Watanabe, Phys. Rev. A 81 , 033415 (2010). 

[33] M. B. Gaarde, P. Antoine, A. L’Huillier, K. J. Schafer, and K. C. Kulander, Phys. Rev. A 
57 , 4553 (1998). 

[34] J. Tate, T. Auguste, H. G. Muller, P. Salieres, P. Agostini, and L. F. DiMauro, Phys. Rev. 
Lett. 98 , 013901 (2007). 

[35] M. F. Ciappina, J. Biegert, R. Quidant, and M. Lewenstein, Phys. Rev. A 85 , 033828 (2012). 

[36] K. C. Kulander, Phys. Rev. A 35 , 445 (1987). 

[37] H. G. Muller, Phys. Rev. A 60, 1341 (1999). 

[38] D. Bauer, Phys. Rev. Lett. 94, 113001 (2005). 


27 



[39] C. 1. Blaga, F. Catoire, P. Colosimo, G. G. Paulus, H. G. Muller, P. Agostini, and L. F. 
DiMauro, Nature Phys. 5 , 335 (2009). 

[40] M. Ruf, H. Bauke, and C. H. Keitel, J. Comput. Phys. 228, 9092 (2009). 

[41] L. V. Keldysh, Sov. Phys. JETP 20, 1307 (1965). 

[42] M. V. Ammosov, N. B. Delone, and V. P. Krainov, Sov. Phys. JETP 64, 1191 (1986). 

[43] J. Grochmalicki, J. R. Kuklinski, and M. Lewenstein, J. Phys. B 19 , 3649 (1986). 

[44] E. Ehlotzky, Nuovo Cimento. 14, 517 (1992). 

[45] M. Lewenstein, P. Balcou, M. Y. Ivanov, A. L’Huillier, and P. B. Corkum, Phys. Rev. A 49 , 
2117 (1994). 

[46] M. Y. Kuchiev and V. N. Ostrovsky, Phys. Rev. A 60, 3111 (1999). 

[47] P. B. Corkum, Phys. Rev. Lett. 71 , 1994 (1993). 

[48] J. L. Krause, K. J. Schafer, and K. C. Kulander, Phys. Rev. Lett. 68, 3535 (1992). 

[49] K. C. Kulander, K. J. Schafer, and K. L. Krause, Super-Intense Laser Atom Physies^ edited 
by B. Piraux, A. L’Huillier, and K. Rzazewski, NATO Advanced Studies Institute Series B: 
Physics, Vol. 316 (Plenum, New York, 1993) p. 95. 

[50] O. Smirnova and M. Ivanov, Attoseeond and XUV Physies: Ultrafast Dynamies and Spee- 
troseopy, edited by T. Schultz and M. Vrakking (Wiley-VCH Verlag GmbH &; Co. KGaA, 
Weinheim, Germany, 2014). 

[51] M. Y. Kuchiev, Sov. Phys. JETP Lett. 45 , 404 (1987). 

[52] F. Brunei, Phys. Rev. Lett. 59 , 52 (1987). 

[53] F. Brunei, J. Opt. Soc. Am. 7, 521 (1990). 

[54] P. B. Corkum, N. H. Burnett, and F. Brunei, Phys. Rev. Lett. 62, 1259 (1989). 

[55] L. Elsgoltz, Differential Equations and the Caleulus of Variations (University Press of the 
Pacific, Miami, United States, 2003). 

[56] H. Goldstein, C. Poole, and J. Safko, Classieal Meehanies^ 3rd edn. (Addison Wesley, San 
Francisco, United States, 2002) p. 356. 

[57] P. B. Corkum, Phys. Rev. Lett. 71 (1994). 

[58] B. Chapman, G. Jost, and R. V. D. Pas, Using OpenMP:Portable Shared Memory Parallel 
Programming (The MIT Press, Cambridge, MA, United States, 2007). 

[59] M. D. Feit, J. A. Fleck, and A. Steiger, J. Comput. Phys. 47, 412 (1982). 

[60] M. Frigo and S. G. Johnson, “Fftw,” http://www.fftw.org/. 


28 


[61] D. G. Arbo, S. Yoshida, E. Persson, K. 1. Dimitriou, and J. Burgdorfer, Phys. Rev. Lett. 96, 
143003 (2006). 

[62] M. F. Kling, J. Rauschenberger, A. J. Verhoef, E. Hasovic, T. Uphues, D. B. Milosevic, H. G. 
Muller, and M. J. J. Vrakking, New J. Phys. 10, 025024 (2008). 

[63] T. Wittmann, B. Horvath, W. Helml, M. G. Schatzel, X. Gu, A. L. Cavalieri, G. G. Panins, 
and R. Kienberger, Nature Phys. 5, 357 (2009). 


29 



